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which the spin-flop transition (a first-order quantum phase transition) 
happens in the antiferromagnetic environment, the decoherence of the 
central spin reaches its highest point. This result is consistent with sev- 
eral recent quantum phase transition witness studies. The influences 
of the environmental temperature on the decoherence behavior of the 
central spin are also investigated. 
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I. Introduction 



Maintaining sufficient quantum coherence and quantum superposition 
properties is one of the most important requirements for applications in 
quantum information processing, such as quantum computing, quantum tele- 
portation [1,2], quantum cryptography [3], quantum dense coding [4], and 
telecloning [5]. Decoherence is a process by which a quantum superposition 
state decays into a classical, statistical mixture of state. It arises from the un- 
avoidable interaction and thus entanglement between quantum systems and 
the environments. When the environmental degrees of freedom are traced 
out or averaged over, the quantum states of the quantum systems are no 
longer pure and the quantum systems are said to have decohered. 

In recent years, the decoherence behavior of a small quantum system, par- 
ticularly a spin system, interacting with a large environment has attracted 
extensive attentions [6-17]. Different environments result in different deco- 
herence features [18]. In some cases of the problems of a central spin coupled 
to an environment, a spin environment may be more appropriate to rep- 
resent the actual localized background spins or magnetic defects than the 
delocalized oscillaor bath that is usually used. The spin environment could 
be lattice nuclear spins, Ising spin bath, spin-star environment or the device 
substrates doped with spin impurities. In Ref . [19] , the decoherence behavior 
of a localized electron spin in the nuclear spin environment is investigated 
experimentally. The influence of an external magnetic field is also considered. 
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The decoherence factor is shown to be a Gaussian type. By attributing the 
electron spin decoherence to quantum fluctuations of lattice nuclear spins, 
Ref. [20] attempts to explain the experimental results. There, the signifi- 
cant coupling between the electron spin and nuclear spins is assumed to be 
the Ising-type interaction and the coupling among nuclear spins is, however, 
taken to be the dipole-dipole interaction. Recently, the dynamics of the re- 
duced density matrix for central spin systems in a spin bath described by 
the transverse Ising model has been analyzed using a perturbative expansion 
method [16] or a mean-field approximation [17,21]. The interaction between 
the central spin system and the spin bath for these case was assumed to be 
of a Ising type. It has also been reported recently that the reduced dynam- 
ics and for one- and two-spin-qubit systems in a spin-star environment [22] 
has been analyzed [23-26]. There, the interaction between the system and 
environment was assumed to be of a Heisenberg XY interaction. 

Antiferromagnetic materials have been reported recently to have appli- 
cations in the area of quantum information processing [27-32] . It was shown 
that quantum computing is possible with a wide variety of clusters assemble 
from antiferromagnetically coupled spins [28]. An antiferromagnetic molec- 
ular ring was also suggested [32] to be a suitable candidate for the qubit 
implementation. The decoherence of the cluster-spin degrees of freedom in 
the antiferromagnetic ring is expected to arise mainly from the hyperfine 
coupling with the nuclear spins. Reference [30] investigated the electron spin 
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decoherence rate of the quantum tunneling of the Neel vector in a ring type 
antiferromagnetic molecular magnet through the study of a nuclear spin cou- 
pled to an electron spin in the ring. In this paper, we study the opposite case. 
We investigate the decoherence behavior of a central spin coupled to an anti- 
ferromagnetic environment. This decoherence problem of a central spin in an 
antiferromagnetic environment has been investigated in Ref. (33]- A similar 
problem of a spin-1/2 impurity embedded in an antiferromagnetic environ- 
ment has been studied in the context of quantum frustration of decoherence 
in Ref. [34J. There the impurity spin is coupled in real space locally to just 
one spin of the antiferromagnetic environment, in contrast to the central spin 
model where the central spin is coupled isotropically to all the spins of the 
environment. In Ref. [35], a large-spin impurity in a ferromagnetic bath is 
investigated in the context of quasi-classical partial frustration of dissipation. 

In the present paper, we consider the case of how an external magnetic 
field will affect the decoherence behavior of a central spin in an antiferromag- 
netic environment. One of the interesting phenomena in antiferromagnetic 
materials under an applied magnetic field is the magnetic-field-induced spin- 
flop transition. When the applied magnetic field is increased to the critical 
field point, the antiferromagnetic polarization flips into the direction perpen- 
dicular to the field. This is called the spin-flop transition. The phenomena 
of the spin-flop transition have been observed experimentally [361 EZj- It is 
thus particularly interesting to investigate how the central spin coherence 
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changes as a function of a globally applied external field in an antiferromag- 
netic environment, especially when the magnetic field strength is near the 
critical field of the flip-flop transition. Indeed, the study of the relationship 
between the pure-mixed state transition of a central spin and the quantum 
phase transition (QPT) of a bath to which the central spin is coupled has 
attracted much attention recently [Ml ESI HD]- A QPT is driven only by 
quantum fluctuations (e.g., happens at zero temperature where the thermal 
fluctuations vanish) and can be triggered by a change of some of the cou- 
pling constants (e.g., external fields) defining the system's Hamitonian. In 
Refs. [38], [391 SO] the concept of Loschmidt echo (or equivalently the deco- 
herence factor) has been used to investigate quantum criticality in the XY 
and anisotropic Heisenberg spin chain models. While the quantum state fi- 
delity approach has also been put forward in Refs. [HI H2] to study QPTs in 
the Dick, BCS-like and XY models. Basically, the signature of QPT (in the 
case of finite spin number) is the dramatic decay of the asymptotic value of 
the Loschmidt echo or the quantum state fidelity at the critical point. That 
is, the closer the bath to the QPT, the smaller the asymptotic value of the 
Loschmidt echo or the decoherence factor of the central spin. The extension 
of the investigations of quantum criticality to finite temperatures has been 
reported in Ref. [12]. We will show in this paper that a similar behavior of 
the decoherence factor of the central spin also appears in our antiferromag- 
netic bath model when the external magnetic field strength is tuned to near 
the critical field point of the flip-flop transition, a first order QPT. 
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It is known that the values of the critical magnetic field of the spin- 
flop transition can be obtained using the spin wave theory [131 0U 05]. We 
thus apply, in this paper, the spin wave approximation to deal with the 
antiferromagnetic environment in a globally applied magnetic field at low 
temperatures and low energy excitations. We obtain explicitly an analytic 
expression of the decoherence factor for the central spin. It is found that 
the external magnetic field affects the decoherence process through its effect 
on the antiferromagnetic environment. Our results explicitly show that the 
decoherence factor which displays a Gaussian decay with time depends on 
the strength of the external magnetic field and the crystal anisotropy field in 
the environment. One of the important results of our present paper is that 
when the magnetic field is increased to the critical field point at which the 
spin-flop transition happens, the decoherence of the central spin reaches its 
highest point, consistent with the QPT witness studies in Ref. [38 | [39l HUl H2] . 
We also investigate the influence of the environmental temperature on the 
decoherence of the central spin. 

The paper is organized as follows. In Sec. II the model Hamiltonian 
is introduced and the spin wave approximation is applied to map the spin 
operators of the antiferromagnetic environment to bosonic operators. In 
Sec. Ill, by tracing over the environmental degrees of freedom, we calculate 
the time evolution of the off-diagonal density matrix elements of the central 
spin and obtain explicitly the decoherence time of the central spin. In Sec. IV, 
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the influence of the environmental temperature, the environmental structure, 
and the external magnetic field on the decoherence behavior of the central 
spin are discussed. Finally a short conclusion is given in Sec. V. 



II. Model and Spin Wave Approximation 

We consider a single central spin coupled to an antiferromagnetic environ- 
ment. It is assumed that the central spin is a spin-| atom and the antiferro- 
magnetic environment consist of atoms with spin S. A global magnetic field 
is applied to both the central spin and the antiferromagnetic environment. 
The total Hamiltonian can be written as 

H = H S + H SB + H B , (1) 

where Hs, H B are the Hamiltonians of the central spin and the environment 
respectively, and H S b is the coupling term [TTJ, [TS], [33J ESI- They can be 
written as 

H s = -gfi B BS z Q , (2) 

HsB = _ ^oE(^« + ^,i)' ( 3 ) 



i,5 j,8 

- 9f M B (B + B A ) ]T SI, - gf i B (B - B A ) £ Sfo (4) 

* 3 

where g is the gyromagnetic factor and \i B is the Bohr magneton. For sim- 
plicity, significant interaction (Eq. 3) between the central spin and the en- 
vironment is assumed to be of the Ising type with J being the coupling 



constant [20]. J is the exchange interaction and is positive for antiferromag- 
netic environment. The effects of the next nearest-neighbor interactions are 
neglected, although they may be important in some real antiferromagnets. 
We assume that the spin structure of the environment may be divided into 
two interpenetrating sublattices a and b with the property that all nearest 
neighbors of an atom on a lie on b, and vice versa S a> j (S^-) represents 
the spin operator of the ith (jth) atom on sublattice a (b). Each sublat- 
tice contains N atoms. The indices % and j label the N atoms, whereas the 
vectors 5 connect atom % or j with its nearest neighbors. B represents a uni- 
form external magnetic field applied in the z direction. The anisotropy field 
Ba is assumed to be positive, which approximates the effect of the crystal 
anisotropy energy, with the property of tending for positive magnetic mo- 
ment fiB to align the spins on sublattice a in the positive z direction and the 
spins on sublattice b in the negative z direction. Reference [10] has analyzed 
the crossover from the case of a single link of the system spin to a bath spin, 
to the case of the central spin model in which the system spin is uniformly 
coupled to all the spins of the bath. This may serve to justify the coupling 
Hamiltonian of Eq. (3). 

We use the Holstein-Primakoff transformation, 

Sa,i = \J%S -afaidi, S~ ti = af^2S - afa { , S* ti = S - afa { , (5) 
S+ = b+y/2S - bpj, S h j x 2S b-bjb,. S^^b+bj-S, (6) 

to map spin operators of the environment onto bosonic operators. We will 
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consider the situation that the environment is in the low-temperature and 
low-excitation limit such that the spin operators in Eqs. (5) and (6) can be 
approximated as ~ \^2Sai, and ~ b^y2S. This can be justified as 
in this limit, the number of excitation is small, and the thermal average < 
afcii > and < bfbi > is expected to be of the order 0(1/ N) and can be safely 
neglected with respected to 25* when N is very large. The Hamiltonians Hsb 
and Hb can then be written in the spin-wave approximation as [HI HU US] 

Hsb = -^S^ibth-ata,), (7) 



H B = E + 2MSJ \ £ af ai + ]T bp, + 2M S J J2( a A+s + a t h ~tJ 

•7 / U 



+gfi B (B + B A ) af^i - 9^b(B - B A ) £ bjbj, (8) 

» 3 

where M is the number of nearest neighbors of an atom and 

E = -2NMS 2 J - 2NSg/i B B A . (9) 

We note here that in obtaining Hamiltonian (8) in line with the approxi- 
mations of S^i ~ y/2Sai, and S^j ~ b^V^S in the low excitation limit, we 
have neglected any term containing products of four operators. The low ex- 
citations correspond to low temperatures, T <C T/v, where Tjy is the Neel 
temperature 



Transforming Eqs. (7) and (8) to the momentum space, we have 

Jo 



H 



SB 



-7=^oE(^-4«k), (io) 



H B = E + (2MSJ + gpL B B A + gpL B B)J2aUk 

k 

+ (2MSJ + gf i B B A - 9l x B B) £ 6+6 k 

k 

+2MSJ5> k (a+fc+ + a k 6 k ), (11) 

k 

where 7 k = M _1 I] ( 5>e lk ' <5 . Then by using the Bogoliubov transformation, 

«k = Mk^k - v k b£ , (12) 
/^k = Mk&k - ^ka k , (13) 



where u£ = (1 + A)/2, ^ = -(1 - A)/2, and A = 1/^1 -7 k , the Hamilto- 
nians -f/ss and H B can be diagonalized as (h — 1) 

^ = -^=^E(^k-<«k), (14) 
ViV 



= ^ + E4 +) (« k V + ^) + E^(^ k + +^ (is) 

where a k (a k ) and /3 k (/3 k ) are the creation (annihilation) operators of the 
two different magnons with wavevector k and frequency co^ (^ k _ ' ) ) respec- 
tively. We note here that the coupling between different magnon modes under 
the spin- wave approximation is neglected. This is a consequence of neglecting 
terms containing products of four operators in Eq. (8) in the low tempera- 
ture limit (T <C T/v) of the spin- wave approximation, where the number of 
excitations of the antiferromagnetic environment is sufficiently small. For a 
cubic crystal system in the small k approximation, 
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where I is the side length of cubic primitive cell of sublattice. E' is a new 
constant. Let k = and ou^ ^ = 0, a critical magnetic field, 



is obtained [I3J HH US] • When the external magnetic field exceeds B c , we have 
a;^ < 0. This indicates that this branch of magnon is no longer stable due 
to the externally applied magnetic field. As a result, the antiferromagnetic 
polarization flips perpendicular to the field, i.e., the magnetic field induced 
spin-flop transition happens. The spin-flop transition demonstrates a signif- 
icant change of the spin configuration in the antiferromagnetic environment. 
This phenomenon has been observed and investigated for many different ma- 
terials 08]. The spin wave theory is known to describe well the 
low-excitation and low-temperature properties of antiferromagnetic materi- 
als. Despite this low-excitation approximation, the spin wave theory also 
describes well the physics for B < B c and the value of the critical magnetic 
field of the spin-flop transition in antiferromagnetic materials [43l 1441 [45] . 
We will thus use the spin wave theory to discuss the decoherence time of the 
central spin under the influence of the antiferromagnetic environment when 
the external magnetic field is tuned to approach B c from below (i.e. B — > B c 
for B < B c ) . It will be shown that an analytic expression for the decoherence 
time can be obtained. 




(17) 



III. Decoherence Time Calculation 
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In this section, we calculate the time evolution of the off-diagonal elements 
of the reduced density matrix for the central spin. We assume that the initial 
density matrix of the total systems is separable, i.e., p(0) = ps(0) <8> ps- The 
initial state of the central spin is described by ps(0). The density matrix 
of the environment satisfies a thermal distribution, that is ps = e~ HB ^ T /Z, 
where Z is the partition function and the Boltzmann constant has been set 
to one. We are interested in the dynamics of the off-diagonal elements of 
the reduced density matrix which are responsible for the coherence of the 
system. This is equivalent to calculating the time evolution of the spin-flip 
operator S = |0)(1|, where |0) and |1) are respectively the lower and upper 
eigenstates of Sq. By tracing out the environmental degrees of freedom, the 
time evolution of Sq = |0)(1| can be written as: 



where the subscript B of tr B indicates that the bath degrees of freedom 
are traced over. We describe briefly how Eq. (18) is obtained. Since the 
Hamiltonian Hb commutes with the total Hamiltonian H, we may move the 
last term in the second line of Eq. (18) to the place next to the operator |0) (1 1 
and evaluate the resultant expression e~* m |0)(l|e* m first using Eqs. (1), (2), 
(14) and (15). After this straightforward evaluation, the operator |0)(1| can 
then be put outside of the trace trg, as shown in the last line of Eq. (18). 



S (t) = trB{e- iHt [S^(0)®p B ]e im } 




|0)(l||tr B [e 



f k -«+a k )t -H B /T 



(18) 
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This is because tr B does not influence the degrees of freedom of the central 
spin. The partition function can be evaluated as follows. 



tr B (e- H */ T ) 



tr, 



"(+)«, 11 , _,.,(-), 



k 



n 



i 



k k k 1 - e-rVT Y i _ g-^-VT 

Similarly, the last part of the numerator of Eq. (18) can be evaluated as 



.(19) 



tr B 



-i glXB Bt-i^L £ k (/3+/3 k -a+a k )t -H B /T 



II- II 

k 



e -ig» B Bt e -E' -Q e -^ +) /2T 
k 



Defining decoherence factor r(t) as 



5 -(0=5 -(0)r(t), 



(21) 



we then obtain from Eqs. (18)-(20), 



where 



r(t) 



y y r -ign B Bt 
+ — ' 
yoyo 



(22) 



yt 



y + 
y~ 



n 



k l-e"4 +) /T' 
?i-e-4-VT' 



n 



i 



= n 



i 



IT 1 - e-^'e-^Vr 
13 



(23) 
(24) 
(25) 
(26) 



To further evaluate the decoherence factor in the thermodynamic limit, we 
proceed as follows. 



lny + = -£ln(l-e-4 +) /^ 
k v ' 

= / ^ (l - e^ +) / T ) Airk 2 dk, (27) 

where V is the volume of the environment. At low temperature such that 
(^k^)max » T, we may extend the upper limit of the integration to infinity. 
With x = kl and iV = V/l 3 , we obtain 

ln?/o + = -^/"^(l-e-^)^. (28) 
In the same way, we have 

l n Vo = -^/ o °°ln(l-e-^)^, (29) 
lny + = --?Lj™ln(l-e ie e-^ +) / T ^x 2 dx, (30) 
lny- = -^j"\n(l-e- i6 e-^~ ) / T ^x 2 dx, (31) 



where 



= (32) 



Letting 



/±(0) = J^ln^l-e^e-^^x'dx 

- j^\n{l-e-^ ±) ' T ^x 2 dx ) (33) 
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we can rewrite the decoherence factor as 



r(t) = e L ^ \^ (34) 



Its absolute value is then 



-Re 



J 2 t 2 



|r(t)| = e L e 2 (35) 

In the thermodynamic limit, i.e., N — > 00, we have, from Eq. (32), — > 0. 
It is then obvious from Eq. (33) that f ± {0) — > as 6 1 — > 0. To find the relation 
between 9 and Re/ ± (6 I ) in the thermodynamic limit, we calculate 

Re/±(0) 



Tj = lim 



e 2 

Relink 

roc g=ie ±ze e k / 2J„ 

JO (±)/T 

Re lim. 1-^-^ 7 



20 

x 2 dx. (36) 



1 roo e-^'/T 



..2. 

(±) . 



2 Jo (1 _ e -^f7T)2 
From the second to the third line of Eq. (136|) . we have used the rule of 
lim I _ lc [/(i)/g(a;)] = lim a .^ c [/'(x)/(7'(a;)], where the prime denotes the deriva- 
tive. The absolute value of the decoherence factor from Eq. (35) in the 
thermodynamic limit (N — > 00, i.e., 6^0) can then be expressed as 

\r(t)\ = e * = e *o, (37) 



where the decoherence time 



to = ^ ■ (38) 
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Equations (37) (38) and (36) are the central result in the present paper. It 
indicates that the decoherence factor displays a Gaussian decay with time. 
The factor t 2 in the exponent is different from the Markovian approximation 
which usually shows a linear decay in time in the exponent. The quantities 
rj in Eq. (36) can in principle be evaluated numerically. An analytic expres- 
sion can, however, be obtained [see Eq. (39)] as the external magnetic field 
approaches the critical field value from the below, i.e., B — > B c and B < B c . 

IV. Results and Discussions 

We discuss the numerical results and figures calculated using Eqs. (36)- 
(38), and present the analytical result when B — > B c from below in this 
section. For a typical antiferromagnetic anisotropy crystal, one has B^ m 
O.lTesla, 2MSJ/gfiB ~ lOOTesla [43J. Our results obtained using the spin 
wave approximation are valid only for not too large temperatures, much 
smaller than the Neel temperature, T/y. According to neutron diffraction 
studies, the Neel temperature of antiferromagnet TbAuIn is 35 K [49] . Some 
antiferromagnets may have higher Neel temperatures. So in the follow- 
ing analysis, the environmental temperature is restricted below T/gfis = 
2.5Tesla, i.e., T « 3.4 K. 

Figure 1 presents the decoherence time as a function of the external 
magnetic field for different temperatures with a crystal anisotropy field of 
Ba = O.lTesla. We can see that the decoherence time decreases with the 
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increase of the external magnetic field. Figure 2 is similar to Fig. 1, except 
B A = 0.15 Tesla. Comparing Fig. 1 with Fig. 2, we find that the larger crystal 
anisotropy field suppresses the decoherence of the central spin. This field- 
dependent decoherence behavior may be inferred from the effective Hamil- 
tonian, Eqs. (14)-(16). From the interaction Hamiltonian, Eq. (14), we see 
that the larger the difference in the magnon excitation number between the 
magnon /3 k and magnon a k , the stronger the effect of the environment on 
the central spin. When the external magnetic field is zero or absence, the 
two magnons, from Eq. (16), have the same frequency for a given wavevector 
k. At a given temperature, the average thermal excitation number may be 
the same for the two magnons, but the fluctuation in the excitations for each 
individual magnon may not be the same at the same time. As a result, trac- 
ing over the magnon states still causes decoherence on the central spin. If 
the external magnetic field is increased, the magnon frequency c<j k ' decreases 
but uj^ increases. Consequently, the magnon mode /3 k is easier to be excited 
than the magnon mode a k at a given anisotropy field and temperature. This 
then results in a larger magnon excitation number difference and fluctuation, 
and thus a stronger decoherence effect. So the decoherence time decreases 
with the increase of the external magnetic field. The effect of the anisotropy 
field on the decoherence time could be understood in a similar reasoning. 
At a fixed external magnetic field, both the magnon frequencies w k and 
uj^ increase with increasing anisotropy fields. Thus the average excitation 
numbers reduce, so is the fluctuation in their difference. As a consequence, 
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Figure 1: The decoherence time as a function of the external magnetic 
field for different temperatures: T/gfis = 0.8Tesla (solid curve), T/gfiB = 
l.OTesla (dashed curve), and T/g^B = 1.2Tesla (dot-dashed curve). Other 
parameters are M — 6, MJ/gfis = lOOTesla, Ba = O.lOTesla. 

the decoherence time increases with the increasing anisotropy field. 

An alternative way to understand the field-dependent decoherence time 
may be in terms of quantum correlations. There is a kind of tradeoff be- 
tween the external magnetic field and the anisotropy field. The anisotropy 
field makes the antiferromagnets stable. On the other hand, the external 
magnetic field tends to reduce the antiferromagnetic order of the environ- 
ment. Therefore the stronger the external magnetic field, the smaller the 
antiferromagnetic order. On the contrary, the larger the anisotropy field, 
the stronger the correlation of the antiferromagnetic environment. If the 
constituents (spins) of the environment maintain appreciable correlations or 
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Figure 2: Same as Fig. 1, except Ba = 0.15Tesla. 

entanglement between themselves, then there is a restriction on the entan- 
glement between the central spin and the environment [501 EE]- As a con- 
sequence, this sets a restriction on the amount that the central spin may 
decohere [Ml [501 EH E2] • Thus as far as the decoherence of the central spin 
is concerned, the anisotropy field has a similar effect to the exchange in- 
teraction strength between the constituents (spins) of the antiferromagnetic 
environment. Strong intra-environmental interaction results in a strong an- 
tiferromagnetic correlation, thus an effective decoupling of the central spin 
from the environment and a suppression of decoherence [501 El] • Therefore 
the decoherence time increases with the increase of the anisotropy field but 
decreases with the increase of the strength of the external magnetic field. In 
summary, the results shown in Figs. 1 and 2 confirm that strong correlations 
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within the environment suppress the decoherence effects [TB I 1501 loTj 152"]. 

In Ref. [34], a similar problem of a spin-1/2 impurity coupled locally 
(in contrast to uniformly in the central spin model considered here) in real 
space to an antiferromagnetic environment has been studied in the context 
of quantum frustration of decoherence. There, the Hamiltonian is reduced 
to the case of a localized impurity spin coupled with two non-commuting 
spin component operators respectively to two bosonic baths. The quantum 
frustration is referred, in Ref. [34], to the lack of a preferred basis for the 
impurity spin due to the competition of the non-commutativity of the cou- 
pled spin operators to the two baths. As a result, the entanglement of the 
impurity spin with each one of the baths is suppressed by the other, and the 
decoherence phenomenon is thus frustrated. On the other hand, the effective 
Hamiltonian of Eqs. (2), (14) and (15) under the spin-wave approximation 
in our central spin model indicates clearly that the Sq eigenstates are the 
preferred basis of the central spin. It also indicates that the influence of the 
two magnon environments seems to partially cancel each other and lead to 
a situation of less decoherence than a central spin coupled to a single bath. 
Though this may also be regarded as some kind of decoherence frustration, 
its cause in this spirit seems different from the quantum frustration of deco- 
herence discussed in Ref. [34] . 

When the strength of the magnetic field is increased further toward the 
critical point where the spin-flop transition occurs, (B = B c « 4.5Tesla for 
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the parameters used in Fig. 1 and B c pa 5.5Tesla in Fig. 2, the decoherence 
time approaches zero quickly. One can obtain analytically from Eqs. (16), 
(36), and (38) for the decoherence time near the critical field as 

AJn s T a \gu B (B c - B)]* . 

lim J r = ^—M*Ji iy ^ BK 39 

b-*b c (b<b c ) T (gti B B e )i 

We note that this result using spin wave theory is valid for B < B c , i.e., before 

the spin-flop transition. Equation (39) indicates that when B approaches 

B c , the environment exerts a great influence on the central spin so that 

the coherence is destroyed thoroughly. This is a catastrophe for quantum 

computing based on such spin systems at the critical field B c . Our result 

in the thermodynamic limit is consistent with those in Refs. [381 1391 HO, H2] 

for the case of finite spin number N. There, the asymptotic value of the 

Loschmidt echo of a bath (or equivalently the decoherence factor of a central 

spin) is dramatically reduced near the critical point and may serve as a good 

witness of QPT in the case of finite N. In our case, the decoherence of 

the central spin is enhanced near the critical point of the antiferromagnetic 

bath. The closer the bath to the QPT (the spin-flop transition from below, 

B < B c ), the smaller the decoherence time of the central spin (see Figs. 1 and 

2). The decoherence reaches its maximum when the values of the external 

magnetic field is increased to the critical field point of the spin-flop transition. 

The Gaussian decay behavior of the decoherence of a central spin similar to 

Eq. (37) has also been reported and discussed in Refs. [381 1391 HO] . 

We discuss the influence of the environmental temperature on the de- 
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coherence behavior of the central spin below. The decoherence behavior is 
sensitive to the environmental temperature for weak external magnetic fields. 
High temperatures cause strong decoherence (see Fig. 1). When the external 
field is strong, the influence of the environmental temperature becomes minor 
or weak. Finally, the decoherence time is less sensitive to the change of the 
external magnetic field at high temperatures. This can be seen from Fig. 1 
that the low temperature case (solid curve) has a sharper curve dependence 
than the high temperature cases (dash and dot-dash curves). 
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Figure 3: Time evolution of the decoherence factor for different tempera- 
tures: T/gfis = 0.8Tesla (solid curve), T/qhb = l.OTesla (dashed curve), 
and T/g/jbB = 1.2Tesla (dot-dashed curve). Other parameters are M = 6, 
MJ/gfi B = lOOTesla, B = 1 Tesla, B A = O.lOTesla. 

In Fig. 3, we plot the time evolutions of the decoherence factor for dif- 
ferent temperatures with the external magnetic field fixed at B = 1 Tesla. 
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Figure 4: Same as Fig. 3, except B = 2 Tesla. 

As expected, the decoherence factor, under the influence of the environment, 
decreases with time. At a higher temperature, the decay rate is larger as 
shown in the figure. Figure 4 is similar to Fig. 3, except that the external 
magnetic field is at B = 2 Tesla. We see again that the coherent behavior 
of the central spin is suppressed with the increase of the external magnetic 
field. 

V. Conclusion 

We have studied the decoherence of a central spin coupled to an antifer- 
romagnetic environment in the presence of an external magnetic field. The 
results, obtained using the spin wave approximation in the thermodynamic 
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limit, show that the decoherence factor displays a Gaussian decay with time. 
The decoherence time decreases, as expected, with the increase of tempera- 
tures. Furthermore, the external magnetic field promotes decoherence effects. 
The decoherence reaches its highest point at the critical field of the spin-flop 
transition, consistent with the QPT witness studies in Ref . [38l [39j EES H2] ■ 
In contrast, the strong anisotropy field suppress decoherence of the antifer- 
romagnetic environment on the central spin. Therefore, in order to reduce 
the loss of coherence of the the central spin, we could decrease the envi- 
ronmental temperature, eliminate the external magnetic field, and choose 
the antiferromagnetic surrounding or underlying antiferromagnetic materials 
with a strong crystal anisotropy field. 
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